Designed nonlocal pseudopotentials for enhanced transferability 
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A new pseudopotential generation metliod is presented wliicli significantly improves transferabil- 
ity. The method exploits the flexibility contained in the separable Kleinman-Bylander form of the 
nonlocal pseudopotential [Phys. Rev. Lett. 48, 1425 (1982)]. By adjusting the functional form 
of the local potential, we are able to improve the agreement with all-electron calculations. Results 
are presented for the Ca atomic pseudopotential. Configuration testing, logarithmic derivatives and 
chemical hardness all confirm the accuracy of these new pseudopotentials. 



The pseudopotential (PS) approximation, or the sep- 
aration of electrons into core and valence based on their 
level of participation in chemical bonding, lies at the 
heart of most modern electijonic structure calculations. 
The original PS formalisnxl |grew out of the orthogo- 
nalized plane-wave approach. □ The atomic PS replaces 
the nuclear Coulomb potential plus core electrons, thus 
simplifying the original system of differential equations. 
Adopting the PS approximation may introduce some un- 
physical results if the PSs are not constructed judiciously. 
The accuracy of the PS, or its transferability, may be 
gauged by its ability to reproduce the results of all- 
electron (AE) calculations in a variety of atomic envi- 
ronments. Most transferability testing has centered on 
configuration testing, characterization of the scattering 
properties using logarithraic derivatives, and, more re- 
cently, chemical hardness. B 

The earliest PSs generated for use in density-functional 
theory calculations replaced the strongly attractive 
Coulombic potential near the origin with a weaker local 
potential, and core electrons were eliminated from the 
calculations .□ In this approach, approximate agreement 
between PS and AE eigenvalues and logarithmic deriva- 
tives was achieved for many elements. However, first row 
non-metals and first-row transition metals could not be 
accurately described by these PSs. 

To improve PS transferability, more complicated (and 
more flexible) semilocal (SL) PSs were designecB with a 
different spherically symmetric potential for each angular 
momentum. This added flexibility permits the enforce- 
ment of the norm-conservation condition at the reference 
energy, e^, for R greater than the core radius, 
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where 4>f-^ (r ) and ^ (r ) are the AE and PS Kohn- 
Sham eigenstates for the state i. Including this criterion 
into PS generation greatly improves transferability. 

Incorporating norm-conservation makes it possible to 
have exact agreement between the AE and PS eigenval- 
ues and wave functions outside for one electronic ref- 
erence state (RS)J3 However, the corresponding single- 
particle differential equation for a PS constructed with 



this method is more complicated because of the angular 
momentum projection. Expressing the SL PSs within a 
plane-wave basis requires the computation of V(G, G') 
instead of just V{G — G') where G is a reciprocal lattice 
vector. This results in a huge memory expense. 

The fuily separable nonlocal (NL) Kleinman-Bylander 
PS formtl dramatically reduces the memory cost of the 
SL PSs. These PSs are constructed from a local po- 
tential and angular-momentum-dependent NL projec- 
tors. In Fourier space, the projector can be expressed 
as W{G) ■ W{G') replacing V{G,G'). This reduces the 
PS memory scaling from iV^ to N. With the inclusion 
of the NL projectors, the resulting single-particle Kohn- 
Sham equation becomes an integrodifferential equation. 
The solutions to an integrodifferential equation may vi- 
olate the Wronskian theorem and possess noded eigen- 
states lower in energy than the nodeless solution.!] A 
simple diagnostic procedureQ allows for detection of these 
lower-energy or ghost states. The separable form of these 
potentials also permits efficient evaluation in solid-stata 
calculations with or A''^ logj CPU-time scalin^'E£l 
for the NL energy contribution and its gradients. These 
potentials have proven very effective for the-study of com- 
putationally intensive large-scale systems.llJ 

While configuration testing and characterization of the 
scattering properties are important tools in ascertaining 
the transferability of a PS, there are some properties that 
cannot be sampled using these diagnostics. These re- 
maining properties involve effects of electrostatic screen- 
ing and non-linearity of the exchange-correlation energy. 
Recently, chemical hardness conservation has been used 
as an effective measure of how these self-consistent terms 
vary with electronic configuration and as an indication 
of accurate PS generation.^ The chemical hardness is de- 
fined as the matrix 
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where Etot[p (t^)] is the total energy of the atom and is a 
functional of the total electronic charge density, p (r ), /„; 
is the occupation number of the n/-th state and e„; is the 
self-consistent nl-th eigenvalue. In the second equality. 
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we have used the fact that e„; = dEtot[p {r)]/dfni- In or- 
der to improve the chemical hardness testing results, the 
form of the local potential has also been studied. Local 
potentials were constructed from iiarious linear combi- 
nations of the semilocal potentials .llHl Chemical hardness 
testing showed that the accuracy of these PSs approached 
the accuracy of a SL PS from the same set of /-dependent 
potentials, but did not exceed it. 

One of the other major objectives in PS generation, be- 
sides transferability, is rapid convergence in a plane-wave 
basis. It has been shown that the residual kinetic energy 
of the RS pscudo-wavefunctions lying beyond the plane- 
wave cutoff energy is an excellent predictor of the basis, 
set convergence error of the PS in a solid or molecule. t2l 
The optimized PS construction is designed to minimize 
this residual kinetic energy. _ 

The ultra-soft PS constructionliJ was formulated to 
generate highly transferable PSs with rapid convergence 
in a plane-wave basis. To facilitate these improvements 
additional NL projectors for additional nl states are in- 
cluded in the form of the PS. The addition of multiple 
projectors per angular momentum allows exact agree- 
ment between the PS logarithmic derivatives and their 
slopes with the AE results for more than one energy. In 
order to improve pseudo-wavefunction smoothness, the 
norm-conservation constraint on the wavefunction is re- 
laxed. To re-introduce norm conservation, a compensat- 
ing valence charge density is added. 

In this paper we present a method for NL PS construc- 
tion according to the KB separable form which improves 
accuracy while retaining the convenience of a single- 
projector representation. 

To construct a PS, an electronic RS is chosen, and 
an AE calculation is performed. Then, a PS, Vps, and 
pseudo-wavefunctions, 0!^f ), are chosen which satisfy 
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where T is the single-particle kinetic energy operator 
and Vn [p] and Vxc [p] are the self-consistent Hartree and 
exchange-correlation energy operators, respectively. The 
latter two operators are functionals of the total charge 

density, p{f), where p{r) = Y.ni fm l^t^nf (^)| • We re- 
quire that the pseudo-wavefunctions obey the following 
criteria: 
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If Vps is not /-dependent, the resulting PS is called local 
(radially and angularly local). More generally, Vps can 
be separated into a local potential and a sum of short- 
ranged corrections: 



where 
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For a SL PS, the corrections, AVf^^, are projection op- 
erators in the angular coordinates and local in the radial 
coordinate. To construct a fully separable NL PS, AVJ is 
formed according to: 



When Vps operates on \(j)^f) we obtain 
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Therefore for the RS we are guaranteed exact agreement 
between the eigenvalues and wavefunctions (beyond 
of the SL and NL atoms. 

However, for any state, other than the RS: 



(8) 



where v! is not required to equal n. The inequality in 
equation (8) illustrates the difhculties involved in assess- 
ing and improving the transferability of NL PSs: the 
transferability of a NL PS can be dramatically different 
from the corresponding SL PS. 

For simplicity, we have focused on the single-projector 
NL PS construction in the current approach. We believe 
that implementation of the presented method along with 
the ultra-soft construction and multiple projectors, may 
provide even greater transferability. We begin by con- 
structing an optimized SL PS. Since equation (4) is sim- 
ply an addition of local and NL terms, we may alter 
arbitrarily without losing the exact agreement between 
the AE and NL eigenvalues and pseudo-wavefunctions at 
the RS, provided we adjust the NL corrections accord- 
ingly. However, we change the eigenvalue agreement at 
any other configuration by doing so. 

Operationally, an additional electronic configuration or 
design state (DS) is chosen. A local augmentation opera- 
tor (A) is added to the local potential forming a designed 
nonlocal (DNL) PS. The augmentation operator is sub- 
tracted from the NL corrections AV; in the following way: 



Vps = (y'°' + 1) + ^ Ai^ 
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TABLE I. Configuration testing for the Ca atom. Eigenvalues and ABtot are given for an all-electron (AE) potential. The 
reference state is Ss'Sp^As^Sd" . Absolute errors are given for a semilocal (SL) pseudopotential, an undesigned nonlocal (UNL) 
pseudopotential and a designed nonlocal (DNL) pseudopotential generated with the presented method. All energies are in Ry. 
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When the augmented Vps operates on the RS, the result 
reduces to equation (7). However for any state other than 
the RS, the second term in equation (9) will contribute 

differently. 

By adjusting A, it is possible to obtain almost exact 
agreement between the AE and DNL eigenvalues for a 
particular DS. With the proper selection of DS electronic 
configuration, excellent transferability can be obtained 
for a variety of ionized and excited states. 

For the Ca atom, we have chosen a RS which is highly 
ionized and a DS which is neutral. All atomic energy cal- 
culations were done within the LDA and optimized PS 
generation methods were used. A 50 Ry plane-wave cut- 
off energy was used to achieve excellent plane- wave con- 
vergence of the pseudo-wavefunctions for the atom. We 
have chosen the s angular momentum channel to be the 
foundation for the local potential and a square step po- 
tential barrier as our augmentation operator. The height 
and radius of the step have been adjusted to reproduce 
the AE eigenvalues for the DS. 

Our selection of RS deserves additional mention. For 
Ca, semi-core states were included as valence in the PS. 
The inclusion of these states allowed for the removal of 
ghost levels, greater local potential design flexibility and 
better overall transferability of the PS. It is also impor- 
tant to note that although we have included multiple s- 
channel states, we do not treat these states with different 
projection operators. The PSs are generated using only 
one NL projector for each angular momentum channel. 

The RS for the Ca PS is 3s^3p^4:S°3dP. The core re- 
gions extend to 1.29, 1.60 and 1.27 a.u. for the s, p and 
d states, respectively. We selected a neutral electronic 



configuration (3s^3p^4s^3d°) as our DS. We found that 
a square potential step with height of 6.76 Ry and radius 
of 0.93 a.u. as the augmentation operator gives excellent 
transferability. 

Table I contains configurational testing results for the 
Ca atom. We present results from an AE atom, a SL PS, 
an undesigned NL (UNL) PS generated with an unaug- 
mented local potential consisting of only the s angular 
momentum channel potential, and a DNL PS which in- 
cludes the square step. The potentials were tested in both 
ionized and excited electronic configurations. In addition 
to eigenvalue information, the table also compares total 
energy differences from the RS. For our choice of augmen- 
tation operator, we find excellent agreement in eigenval- 
ues and total energy differences. In all configurations 
shown, the errors are below 1 mRy for the DNL results. 
The eigenvalue and total energy errors for the DNL po- 
tential are one to two orders of magnitude smaller than 
the SL and UNL results. 

Figure 1 shows the logarithmic derivative differences 
between the AE potential and all three PS generation 
methods for the s, p and d potentials of the Ca atom. For 
the s and p channels, we find excellent agreement over a 
large energy range between the PSs generated with the 
DNL method and the AE potentials. The energy errors 
in the s and p states for the SL and UNL PSs presented 
in Table I are directly related to the logarithmic deriva- 
tive differences. Interestingly, we find that the DNL PS 
d state logarithmic derivative differs from the AE results 
more than the other methods for the RS. This finding 
suggests that hardness testing may be a better predic- 
tor of transferability than logarithmic derivative testing, 
and we plan to investigate this point further. In anal- 



3 



TABLE II. Chemical hardness testing for the Ca atom. Absolute hardness values are compared for an all-electron (AE), 
semilocal pseudopotential (SL), an undesigned nonlocal (UNL) pseudopotential and a designed nonlocal (DNL) pseudopotential 

generated with the presented method. Hardness values were determined for two different electronic configTiratioiis. Each clement 
of the symmetric hardness matrix, Hni^n'V > is indexed by the change in the nZ-th eigenvalue for a change of the n'l'-th occupation 
number. Units are in Ry per occupation number. 
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ogous logarithmic derivative testing at other electronic 
configurations, the DNL PS d state scattering properties 
are more accurate than the other PSs in the range of 
energies shown in Table I. 

Table II contains chemical hardness testing results for 
Ca AE and pseudoatoms, for two electronic configura- 
tions. In both cases, we find excellent agreement in 
chemical hardness between the AE and DNL potentials 
for each electronic configuration. With errors below 0.5 
mRy, the DNL PS is one to two orders of magnitude more 
accurate than the SL and UNL PSs. 

In this paper, we have developed and implemented a 
fully nonlocal pseudopotential approach using the separa- 
ble form of Kleinman and Bylander. We have exploited 
the implicit flexibility contained within the separation 
of the pseudopotential into local and nonlocal parts by 
including an augmentation operator into the local and 
nonlocal parts of the potential. By adjusting the aug- 
mentation operator, we have been able to achieve almost 
exact agreement with all-electron results for a variety of 
ionized and excited states. In addition to configuration 
testing, we have presented logarithmic derivatives and 
chemical hardness tests. All the tests demonstrate signifi- 
cant improvement over semi-local and standard nonlocal 
potentials. Furthermore, we have shown that for elec- 
tronic configurations that contain multiple states with 
the same angular momentum, it is possible to construct 
a pseudopotential with a single nonlocal projector that 
will yield very accurate results. 
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FIG. 1. DifFcrcncos in logarithmic derivatives between the Ca pseudopotcutial and the all-electron calculation for the 
s (top), p (middle) and d (bottom) states at the reference configuration (3s^3p''4s"3d°). Results are given for a semilocal 
pseudopotential (dotted line), an undesigned nonlocal pseudopotential (dashed line) and a designed nonlocal pseudopotential 
generated with the presented method (solid line). 
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